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Abstract. The emission intensity of diamond samples containing nitrogen- 
vacancy centres are measured as a function of magnetic field along a (111) 
direction for various temperatures. At low temperatures the responses are 
sample and stress dependent and can be modeled in terms of the previous 
understanding of the '^E excited state fine structure which is strain dependent. 
At room temperature the responses are largely sample and stress independent, 
and modeling involves invoking a strain independent excited state with a single 
zero field splitting of 1.42 GHz. The change in behaviour is attributed to a 
temperature dependent averaging process over the components of the excited state 
orbital doublet. It decouples orbit and spin and at high temperature the spin levels 
become independent of any orbit splitting. Thus the models can be reconciled and 
the parameters for low and high temperatures are shown to be consistent. 
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1. Introduction 

There is interest in using the NV^ centre in diamond for apphcations in the 
area of quantum information processing [1], quantum computing [21 El HI [5], and 
magnetometry [SI [7]. The apphcations rely on the spin in the ground state, and 
the behavior of this spin in its environment is largely understood [H [H [9l [10] . The 
excited state is also used, for initialization and for readout, and an understanding of 
spin in this state has been obtained from low temperature excitation measurements 
[llj . The state is an orbital doublet, which is usually split by internal strain, and the 
excitation measurements examined the fine structure in both orbital branches. More 
recently, however, fine structure has also been determined from room temperature 
ODMR measurements [121 [1] and it is not clear how these values related to those 
obtained for the two branches at low temperatures. In this work it is shown that the 
room temperature results do not correspond to either branch, but rather are due to 
an averaging of both branches. 

The excited state fine structure has been obtained from high resolution optical 
excitation measurements of single NV centres [11] [13] and from two laser hole burning 
[141 Il5l fT6] . These techniques are applicable at low temperatures but can not be 
extended to higher temperatures. On the other hand, spin information has also been 
obtained from optically detected magnetic resonance (ODMR) of single centres [T2lll7j. 
In this case the observations were made at room temperature and the techniques 
are not readily extended to low temperature. Instead of these techniques, we use 
level crossing of the spin states as indicated by changes of emission as a function of 
magnetic field. This approach is applicable over the entire temperature range from 
room temperature to helium temperatures, and does not requrie the same level of 
specialist equipment. Room temperature responses have been reported previously "18], 
and the close relationship between ensemble and single site results suggested that a 
study of concentrated crystals could be fruitful. Ensemble measurements using crystals 
allow uniaxial stress to be introduced; and so the studies can involve variation of the 
temperature, the magnetic field and the stress. The focus of this investigation was how 
the averaging process changes the observation between high and low temperatures. 

2. Emission Variation v Magnetic field: Room Temperature 

A sample (designated L) with an NV^ concentration of 10~^^ ppm was aligned with 
a (111) face parallel to the magnetic field direction. The sample was excited using 
a 532 nm wavelength laser and the red emission at right angles detected using a Si 
diode. This emission was monitored as a function of magnetic field and is shown in 
Figure [B^a) . There is about a 10% drop in intensity between gauss and 400 gauss 
followed by a flatter response, with sharp features at 510 ±5 gauss and 1028 ±3 gauss. 
The sharp dip at 1028 gauss occurs when there is an avoided crossing of the ground 
state spin levels for centres aligned with the magnetic field (1^ [T5j . This feature 
gets narrower with closer alignment of the (111) axis with the magnetic field [21], and 
adjustment was made to minimize the width. Emission can also be quenched by radio 
waves. The minimum frequency (20 MHz) at which this could be achieved indicated 
that the transverse field was about three gauss, and hence the misalignment was less 
than 0.2 degrees (taken to be degrees). The emission responses as a function of 
magnetic field (magnetic spectra) were then measured for angles of 2, 4 and 6 degrees 
and these are also shown in Figure[lja). The feature of primary interest is the decrease 
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Figure 1. (a) Measured magnetic spectra at room temperature for various 
misalignment angles as indicated in the common legend. The fine structure on 
the 5f0 G feature in the experimental trace for 0° misalignment is due to spin 
cross-relaxation with nearby nitrogen impurities, and is not of interest here, (b) 
Calculated magnetic spectra for the same misalignment angles, simulating an 
ensemble sample with four different NV~ orientations. The variation in emission 
intensity was calculated using rate equations for the 7-level model shown in Figure 
|3] The upper trace shows the effect of spin mixing (spin level avoided-crossing) in 
the ground state at 1028 G, and the other traces also include similar spin mixing 
in the excited state at 510 G. 

at 510 gauss, but it is complicated by cross relaxation with single interstitial nitrogen 
occurring at the same field value (see Section As can be seen in Figure the cross 
relaxation feature is not present when the sample (N) contains few substitutional 
nitrogen impurities. The emission of this sample was low and not convenient for the 
more detailed studies. 

Uniaxial stress of up to 0.5 GP could be applied using a piston arrangement. It 
was apphed along the (Oil) direction and transverse to the centres giving rise to 1028 
gauss and 510 gauss features. The results are shown in Figure [3] and it can be seen 
that the stress has no affect on the room temperature response. 
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Figure 2. Room temperature magnetic spectra for a sample with low 
concentrations of interstitial nitrogen. The fine structure at 510 G is absent even 
for a well-aligned field. 
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Figure 3. Room temperature magnetic spectra (as in Figure [Tl^a)) with and 
without 0.35 CP uniaxial strain applied perpendicular to the (ill) crystal axis 
with which the magnetic field was aligned. 

3. Interpretation of Room Temperature Response 

The optical transition associated with the nitrogen-vacancy defect is known to be due 
to an E transition at a site of trigonal symmetry. Neglecting the orbital 

doublet nature of the excited state, the dynamics have been described in terms of 
a 7 level model; three spin levels in the ground state, three in the excited state, 
and an intermediate singlet level [22] (The TOj = ±1 spin levels in the ground and 
excited states can be degenerate and if these are treated as single levels the system 
can be considered as a 5 level system). These energy levels are shown schematically 
in Figure m with transition rate parameters which have been obtained from lifetime 
measurements and other transient changes in the emission 22J. The dominant terms 
are those associated with the optical transitions where spin projection is conserved, 
and those associated with decay via the singlet. The singlet path results in a fraction of 
the = ±1 excited spins relaxing to the ground = without visible emission, and 
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Figure 4. Energy level models used for calculations, with transition rate 
parameters in units of 10^s~^ from 1221 . Non spin-conserving transitions were 
included in calculations to give realistic spin polarisation levels, but are left off 
this diagram for clarity. The structure in the intermediate singlets is not important 
here since the lifetime of the ^E' level is very short |23) , and so the system enclosed 
in the dashed box can be treated as a single level. The 7-level model includes the 
three ground spin levels, a single intermediate level and three spin excited state 
levels (and is shown in blue). To describe the low temperature results, it is 
necessary to add the three spin levels for the second orbital branch in the excited 
state (shown in orange) to give a 10-level model. 



as a consequence the emission associated with mg = ±1 is weaker than that associated 
with TOs = 0. Once the system is cycled several times the population is predominantly 
in the vris — level (the system is "spin polarized") and the emission is high. Spin 
polarization and emission will be reduced whenever the spin states are mixed, and this 
can occur with the application of an external magnetic field. The mixing introduced 
by a magnetic field is readily calculated and a rate equation treatment of the 7 levels 
can be used to predict the changes in the intensity of the emission as a function of the 
magnetic field (Figure [TJb)). 

The largest changes occur for magnetic fields giving rise to crossing of the spin 
levels. For example, the zero field separation between the TOs = ±1 and TOs = spin 
states in the ground state is 2.88 GHz and so if an axial field of 1028 gauss [g — 2) 
is applied the = and to^ — —1 become degenerate. Any additional transverse 
magnetic field lifts this degeneracy and the basis states become totally mixed. Spin 
polarization is no longer maintained and the emission intensity is reduced. This 
accounts for why there is a sharp reduction in the emission intensity when a magnetic 
field aligned close to a (111) direction is swept in magnitude through 1028 gauss. With 
misalignment, and thus a larger transverse field, the spin mixing occurs over a larger 
range of field and accounts for the broadening of the features. For ensemble samples 
the non-axial centres also contribute to the emission and there is an additional change 
in emission as a function of magnetic field. For a field sweep from to 1400 gauss 
calculation gives a gradual decrease in emission intensity plus the sharp dip at 1028 
gauss. Such emission changes are observed, but the calculation only allowing for this 
ground state avoided crossing does not give good correspondence with the observations 
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as can be seen by comparing upper trace in Figure \V[h) with experiment in Figure 

ma). 

In the present work we also consider the effect of a zero field splitting in the excited 
state. To account for a level crossing at 510 ± 5 G, the zero field splitting is taken to 
be 1.43 ± 0.01 GHz (with g = 2). The effect on the eigenstates can be calculated in 
the same ways as for the ground state. The excited state avoided crossing causes a 
decrease at the field value of 510 gauss which again broadens with misalignment. Also 
the contribution for the non-axial centres gives reduced emission at low field values. 
With this extended rate equation model the combined effect of the ground and excited 
state avoided crossings is found to give excellent agreement with experiment (Figure 

m- 

It is concluded that the room temperature magnetic field data can be explained 
assuming there is an excited state with a single avoided crossing at 510 gauss. Within 
experimental uncertainty, this is consistent with an state having a zero field splitting 
of 1.42 GHz and g = 2 as has been observed in the room temperature ODMR 
measurements [12l [5]. 

4. Other Features 

The interpretation for well aligned fields must be considered more carefully as the rate 
equation model used above is unreliable. In the model it has been assumed that the 
optical transitions occur between the eigenstates of the ground and excited state but 
time scales can become too short for the system to relax to its eigenstates. Consider, 
for example, the situation immediately after there is decay from the singlet into the 
rUg = level in the ground state. Where there is a ground state avoided crossing, 
the effective magnetic fields are all weak (only the weak transverse field is effective) 
and so the spin projection precesses slowly. With intense excitation, the system can 
be pumped out of the ground state before there is any reorientation of the spin. The 
spin polarization and high emission will be maintained. This is inconsistent with 
excitation out of the eigenstates which would conclude that rng = and = — 1 are 
equally populated giving no spin polarization and low emission. Similarly if there is 
avoided crossing in the excited state the spontaneous decay can occur before there is 
reorientation, and again emission remains high inconsistent with the result of emission 
from the eigenstates. Whenever the misalignements are large the reorientation of the 
spin projections occurs quickly compared to the lifetimes, and the time average values 
gives the same result as with the eigenstate description. The calculations used above 
are valid for misalignments greater than 1 degree, and this is the range considered by 
the experiments used in this paper. 

As just discussed, when the magnetic field is aligned within 1 degree the axial 
centre maintains good spin polarization even at 510 gauss. However, at this field value 
the spin polarization can be quenched due to a cross relaxation process. For example 
at room temperature at 510 gauss there is an energy match of the separation of the 
spins of the NV~ centre and S = ^ centres such as substitutional nitrogens. Due to 
the energy match the optically induced spin polarization can not then be maintained 
in the presence of spins with reversed spin population. This cross relaxation including 
the involvement of changes of nuclear spin projection has been the subject of several 
studies [19l[20l[2l[25l[26l[27l[28|. Such a feature can be seen in the 0° misalignment 
trace on Figure [TJa), and it clearly has structure. The feature is observed whenever 
there is a level of spin polarization. This cross relaxation feature is not the subject of 
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Figure 5. Temperature dependance of the magnetic field spectrum. The spectra 
have been fitted beyond f250 G with exponential functions, and normalised at 
the high field limit. The fine structure at 514 G is not related to the excited state, 
but provides useful confirmation of the field magnitude. The inset shows the zero 
field to high field ratio as a function of temperature. 



the present study and adds a level of complexity for the interpretation of the spectra. 

5. Emission Variation v Magnetic field:Low Temperature 

Magnetic field spectra were recorded for temperatures between 8 and 300 K, and 
the variation for sample L is shown in Figure [5j Due to thermal expansion of the 
components holding the diamond sample and thermal variation in absorbtion it was 
difficult to keep the absolute magnitude calibrated over a large temperature range. 
However the emission is determined by the optical pumping cycle, and at high field will 
have a fixed value as the spin eigenstates are (to a good approximation) determined by 
the magnetic field and thus are independent of temperature. The emission responses 
have, therefore, been normalised to their high field intensities. An exponential fit to 
the high field tail of the measured spectra allowed each trace to be extrapolated to its 
asymptotic value. 

The 1028 G feature that corresponds to an avoided crossing in the ground state 
is approximately constant, indicating that the ground state properties do not change 
with temperature. However, the 510 G feature attributed to the excited state changes 
dramatically. At room temperature it is centred on 510 G, slightly narrower than the 
ground state crossing feature, and retains these properties down to 100 K. Cooling 
further to 40 K causes this feature to diminish significantly. By 25 K there is no 
feature at 510 G, and instead there is a local minimum in the spectrum at about 250 
G which is much broader than the ground state crossing. Cooling below 25 K caused 
a double-bump to appear (Figure [2]) . Also it is noted that as the sample is cooled 
the emission intensity at zero- field varies significantly (inset in Figure [5]) . It has a 
minimum at 25 K and a maximum at 100 K. There is a small reduction in zero-field 
intensity as the temperature increases from 100 K to room temperature. 

Unlike the effect at room temperature, uniaxial stress was found to alter the 
magnetic field spectrum at low temperature. Examples are shown in Figure[6] Sample 
L exhibits significant changes with stress. The doublet is changed to single broad 




Figure 6. Low temperature (15 K) magnetic spectra witii and without 0.35 GP 
uniaxial stress perpendicular to the (111) crystal axis with which the magnetic 
field was aligned. The traces for Sample H (higher concentration) have been 
translated up for clarity and correspond to the right-hand vertical axis. 



feature (plus the cross relaxation feature). The second sample H contained a higher 
concentration of NV^ defects and, although barely changed by external stress, already 
had the characteristics of the first sample with stress applied. 

6. Interpretation of Low Temperature Response 

At low temperatures fine structure attributed to the '^E state has been observed in 
holeburning [TH [551 [13 [IS] , excitation [IT1[T3] and photon echo [301 [SI] experiments. 
The Hamiltonian associated with this structure is given by 

H — Hso + 'Hss + "Hstr + 'Hz (1) 

where l-Lso is spin-orbit, Jisa spin-spin, Jistr strain and "Hz the Zeeman interaction. 
The diagonal component of spin-orbit Jiso = A^L^Sz splits the state into three equally 
separated doublets. One doublet (E irreducible representation) is not displaced and 
associated with the Sz states (mg — 0). The state where spin and orbital projections 
are parallel (also E irreducible representation) is displaced down in energy by Az 
and the antiparallel alignments (Ai and A2 states) are displaced up in energy by the 
same magnitude. The off-diagonal spin-orbit A^j^j, (L^S^, -I- Lj,Sj^) has only a minor 
effect as its value is known from previous work to be small (\x,y — 0-2 GHz) [TTj . 
Spin-orbit can mix ^E states with singlets and is important for population decay, but 
has negligible effect on the energy levels and is not included here. Spin-spin Ti^s 
has two contributions. There is a normal spin-spin term Dgs (S^ — |) which shifts 
the vHs = ±1 states with respect to = states by Des- The second spin-spin 
contribution is associated with spins on the two orbital components of the degenerate 
state ^2]. It only affects A\ and Ai states and we assign the magnitude of this 
splitting as A. The energy levels in perfect Ca^ symmetry are therefore shown on 
the left of Figure [7] [TTl [33] . The energies are determined by four parameters and 
the value of the parameters have been obtained from high resolution excitation traces 
of single centres. The most reliable parameters from recent measurements |13) are 
\z = 5.5 GHz, D^s = 1-42 GHz, A = 3.1 GHz, and A^.^ = 0.2 GHz. 
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Figure 7. (a) Excited state energy level splittings as a function of strain, (b) 
Splitting of energy levels for a centre at strain I (0.2 GHz) with a magnetic field 
applied at 2° to the NV~ axis, and (c) the magnetic spectra calculated for this 
strain from the 10-level model in Figure |4] (d) and (e) show the corresponding 
responses to an identical magnetic field when the centre is at strain II (4 GHz). 
The magnetic spectra in (c) and (e) are shown in more detail in Figure [Sfa), but 
are included here to indicate their connection with the spin level avoided crossings. 



The splitting of the excited state energy levels with strain, as shown in Figure 
[7{a), can be calculated by diagonalising the energy matrix 
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where 6x and 5y are strain parameters. 

Strain perpendicular to the trigonal axis lowers the symmetry and Tistr has the 
effect of splitting the into two orbital branches each with three spin levels (Figure 
[TJa)) [TTl[33j- In the lower branch spin-orbit and spin-spin have opposite signs. At 
zero strain spin-orbit is larger resulting in the = ±1 states being lowest. However, 
spin-orbit splitting is reduced with increasing strain and the separation of the m^. = ±1 
and TOs = levels is reduced. The states can cross, although the crossing is changed 
to an avoided crossing by off-diagonal spin-orbit interaction. The distance of closest 
approach of these levels enables the strength of the interaction to be determined and a 
value of Xx,y — 0.2 has been obtained in this way 111]. In the upper branch, spin-orbit 



Time-averaging within the excited state of the nitrogen-vacancy centre in diamond 10 





(a) ' 1 


1 




\ 

\ 0.2 GHz 




0.9 






units) 






9 GHz 




^ 0.8 






sity (ai 




Inten 


(b) ' 1 


1 


a 












E 
W 






0.9 






0.8 








1 , 1 1 1 



500 1000 

Magnetic field (G) 



Figure 8. Calculated magnetic spectra, with a field misalignment of 2°, for (a) 
strain splittings of 0.2, 4 and 9 GHz; and (b) an ensemble witii strain distribution. 
The strain distribution in (b) consisted of 0, 0.2, 2, 4, 7, 9 GHz with corresponding 
weights of 1, 3, 5, 5, 3, 1. The ofT-axis NV~ centres are included in these 
calculations, simulating an ensemble crystal sample. 

and spin-spin have the same sign and the states are always weU separated such that 
there is Httle mixing [Til ES] ■ 

For low fields (< 2500 gauss) the Zeeman interaction is dominated by 2S/3B and 
the splittings are illustrated in Figure [7I^b) and (d) for centres with various values of 
strain. In the case where strain is sufficient to split the ^^E state into two triplets (II) 
each triplet is split analogous to that of an orbital singlet. The important features 
are the avoided crossings of the TOj = and rus = — 1 spin levels. There is one 
associated with each orbital branch. Rate equations for the 10 levels (Figure 2]) are 
used to calculate the changes of emission as a function of field. Decrease in emission 
is obtained at field values corresponding to the two excited state avoided crossings 
plus the one in the ground. As can be seen from the examples magnetic field value 
at which the excited state dips occur varies with strain and can occur at zero field 
(FigureHa)). 

For an ensemble the emission will have contribution from centres with a range 
strain environments, and an example using six strain values is shown in Figure 
[5{b). Reliable calculation would require knowledge of the distribution, but this is 
not available. However, we can comment on how certain characteristics arise. A 
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contribution from centres with mixing at zero field will depress the emission at zero 
field and this has been observed in the experiment (Figure [5]) . Also as the avoided 
crossing in the lower branch always occurs at low field values (< 500 gauss) it is 
understandable why the emission in the 200 - 500 gauss range is low. Conversely 
avoided crossing in the upper branch occur at higher fields and accounts for the slow 
increase in emission intensity as the field is extended above 2000 gauss. 

The conclusion is that the low temperature responses are consistent with that 
expected for the known fine structure of the NV centre and its variation with strain. 
A variation from centre to centre and sample to sample is understandable in this 
model, but this is totally different to the behavior at room temperature. 

7. MCD of Zero-Phonon Line 

An axial magnetic field splits an orbital E state into two components and i?_ 
and these components give allowed transitions to an orbital A state in right and left 
circular polarized light respectively. Measurements of the electronic transition can be 
made using modulation techniques, and by using the first moment of the MCD the 
Zeeman splitting can be determined even in the presence of inhomogeneous broadening 
[34] . The splitting allows the orbital angular momentum to be determined and for 
an A ^ i? would be temperature independent The situation here is slightly more 
complicated as the transition is E rather than d.xi A ^ E. For this case the 

electronic Zeeman term leads to three sets oi A ^ E transitions, one for each spin 
projection [rag = 0,1 and — 1) and the MCD of each has a slightly different magnitude 
[55] . Spin-orbit interaction has the effect of increasing the MCD for vris = — 1 and 
reducing it for rus = +1. The different magnitudes leads to a temperature dependence 
of the MCD. However as spin orbit is only 5.24 GHz the change of magnitude is only 
5 % and much smaller than considered previously |14j . 

The MCD and absorption of the zero-phonon line at 637 nm was measured for 
the 10^^ ppm concentration sample H at various temperatures using a field of 5 Tesla 

The variation in the first moment normalized to the strength of the absorption 
is shown in the insert (Figure [TO)l . The magnitude of the signal indicates a splitting 
of 7 GHz consistent with a g = 0.1 as obtained previously [14]. What is significant 
however, is that the signal drops with increasing temperature and is continuing to 
drop at 200 K. The size of the decrease is much larger the 5 per cent that could be 
attributed to spin-orbit interaction. The interpretation is that the effective angular 
momentum is high at low temperature but approaches zero at room temperature. This 
variation will be discussed below. 

8. Averaging Process With Temperature 

The change of behavior with temperature is attributed to an averaging process and 
will be discussed here. 

It is well known that increasing temperature gives rise to a variation of the 
homogeneous line width of optical transitions. Changes of width involve the transitions 
combining with the absorption and emission of lattice vibrations of slightly different 
energies. The vibrations also induce transfer between electronic states and give the 
standard Boltzman distribution in the population of adjacent states. By considering 
specific optical line widths we are able to estimate the time scale of such processes. 
At helium temperatures a ^A2 E line width of 13 MHz has been observed [36] and 
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Figure 9. Absorption and MCD response for various temperatures. 
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Figure 10. Effective g^s value as a function of temperature, obtained from the 
first moment of the MCD response. 
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is determined by the radiative lifetime of 12 ns. The hfetime does not change with 
temperature and at higher temperature the larger widths must reflect lifetimes of 
the coupled electron-vibrational system. For example at intermediate temperatures 
holeburning gives an estimate of homogeneous line widths and [37] indicate values 
of the order of 20 GHz at 60 K. Above 180 K the transitions are homogeneously 
broadened and can be measured by conventional absorption or emission but give rates 
of 1000 GHz [38 . 

Vibrations interact with electronic orbit but not with spin. Therefore the fast 
rates will involve rapid interchange between orbital components of a degenerate state 
without changing the spin projection. This is what occurs here, and the fast orbital 
change has significant effect on observations. In large magnetic fields the interchange 
will be between states, E^, with positive and negative angular momentum. Rapid 
variation will result in an average state with less or no orbital angular momentum. 
Thus this type of averaging will reduce the effective orbital angular momentum and 
account for the observed reduction of the MCD. 

For strain the interchange will be between orbital component from to Ey with 
the Sz,Sx, Sy spin projections remaining unaltered. Although the projections stay the 
same the energies change (see Figure [7|). Thus normally the ODMR frequency will 
be a few GHz's different for the two orbital branches. However, rapid change of the 
orbital branch at a rate of many lO's GHz will alter the situation and the ODMR 
measurement will give the average frequency rather than two separate frequencies. 
(This is similar to conditions in liquid state NMR.) The averaging effectively quenches 
the contribution of spin-orbit interaction Xz- Also as Sx spin coupled to E^ and Ey 
is distributed between the Ai and A2 states the averaging quenches the effect of this 
splitting A. The remaining term in the Hamiltonian D^s (S^ — |) determines the 
average Sz Sx (and Sz Sy) separation and the ODMR frequency. At high strain 
the spin-spin would include a strain dependent term E^s (S^ — S^) , and the ODMR 
will give two excited state resonances at ± E^s ■ 

The above averaging process results in an 'average' state with a zero field splitting 
determined by the spin-spin interaction D^s (low strain situation). In order to obtain 
a fit for recent low temperature high resolution excitation measurements a zero field 
splitting of 1.42 GHz is required, and this is confirmed from measurements of the 
separation of the averaged energy of Sz,Sx, Sy levels [13] . Thus there is consistency 
between the low and high temperature frequencies, and the averaging accounts for the 
difference between high and low temperature behaviours. 

The averaging is not associated with a static Jahn- Teller effect. For modest 
orbital strain splitting the averaging can result in the appearance of higher symmetry, 
however the splitting varies over a large range from zero to lOOO's GHz. Clearly if 
there is zero strain splitting there is no Jahn- Teller distortion. Also, with very large 
orbital splittings it is observed that the ODMR signal is split and therefore the system 
still has low symmetry. In the case of Jahn- Teller the averaging would result in higher 
symmetry. 

9. Conclusion 

From this work it is concluded that the spin and optical properties of the NV^ centre 
are much simpler at room temperature than has previously been anticipated. Thermal 
vibrations average the orbital components of the excited state and quench all orbital 
contribution to the relative energy of the spin levels. In effect, the excited state 
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has the properties of an orbital singlet. Spin level energies are determined only by 
spin-spin interaction, just as they are in the orbital singlet ground state. This state 
averaging process is most likely a general phenomenon when thermal transition rates 
between orbital branches are large compared to the spin-orbit interaction. The result 
is an effective decoupling of spin and orbit. In the NV^ system, this averaging has a 
prominent effect due to the very small spin-orbit interaction. At low temperature the 
thermal transition rates are slowed and spin-orbit interaction does become significant, 
and so it is necessary to consider both orbital branches of the excited state. 

This averaging has implications for many earlier studies of the NV~ centre. For 
example, it accounts for why previously observed ODMR signals can be explained in 
terms of a spin triplet associated with a non-degenerate electronic state [12j [17] . It 
also justifies our use of a 7-level model both here and in previous publications [22] to 
explain the optical pumping cycle of the NV~ centre at room temperature. 

At room temperature there is little variation in optical properties from site to 
site, which could make it simpler to develop applications. The averaging eliminates 
the need to seek centres with optimum properties, and justifies determining many 
properties from ensemble measurements. Counterintuitively, stronger spin polarization 
is obtained at room temperature than at low temperature, and this is likely to be of 
considerable interest for the development of quantum devices. 
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